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An initially isotropic medium, when subjected to either a magnetic field or a coherent field, 
can induce anisotropy in the medium and can cause the polarization of a probe field to rotate. 
Therefore the rotation of probe polarization, due to magnetic field alone, can be controlled efficiently 
with the use of a coherent control field. We demonstrate this enhancement of the magneto-optical 
rotation (MOR) of a linearly polarized light, by doing detailed calculations on a system with relevant 
transitions j = <-> j = 1 ^ j = 0. 
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An isotropic medium having m-degenerate sub-levels when subjected to a magnetic field exhibits birefringence in 
^ ' its response to a polarized optical field. This is due to the fact that Zeeman splitting of magnetic sub-levels causes 
asymmetry in refractive indices for left and right circular polarization components of the optical field. The result is 
magneto-optical rotation (MOR); i.e., the plane of polarization of the light emerging out of the medium is rotated 
with respect to that of the incident. If x+{x~) represents the susceptibility of the birefringent medium corresponding 
to right (left) circular polarization component of the probe, the rotation angle is given by 

> . 

O. e = nkpl{x--X+), (1) 

(N 

' where kp corresponds to propagation vector of the probe and I is length of the medium along the direction of 
propagation. The susceptibilities x± ^re functions of characteristics of atomic transitions and hence studies of 9 yield 
information on atomic transitions. Extensive literature on MOR exists. Experiments have been carried out with both 
synchrotron and laser sources . Effects of larger probe powers and resulting saturation have also been studied at 
length g^. 

The present study is motivated by the possibility that the susceptibilities x± can be manipulated by the application 
of a control laser [^-pd|]. The control laser modifies the atomic structure and thus by tuning the strength and 
] the frequency of the control laser one can obtain suitable modifications in x±- Wielandy and Gaeta |l2| recently 
■ demonstrated control of polarization state of the probe field in an initially isotropic medium (see also p^Jl^ ). In this 
^ ' paper we study in detail the possibility of the enhancement of MOR produced by the control laser. 

The organization of the paper is the following: in Sec II, we describe the response of an isotropic medium to a 
' linearly polarized weak probe. In Sec. Ill, we calculate x± of the medium in presence of a control field using the density 
• i-H matrix formalism. We show analytically how a circularly polarized control field causes anisotropy in the medium. 

Further in in Sec IV, we present the numerical results - that demonstrate the control laser induced birefringence and 
^ I we show the resulting large enhancement in MOR in the regions identified by the analytical results of Sec. III. We also 
■ - - I show new probe-frequency domains where significant MOR can occur. 

II. RESPONSE OF AN ANISOTROPIC MEDIUM TO A LINEARLY POLARIZED LIGHT 

Let us consider an incident probe field Ei„ with frequency ujp propagating inside an active medium along the 
quantization axis z is 

^™(z, t) = Ep{z)e-'^-'+'^-' + c.c. . (2) 
We resolve the incident field amplitude £p into two circularly polarized a+ and (t_ components 

Sp — c-^Sp-\- + tp—Sp— ^ (3) 



*An article in honor of Prof. M. O. Scully - who led us into many areas of quantum optics and laser physics. 
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where unit polarization vectors e± corresponding to (t± polarization are defined in terms of unit vectors x and y as 

The polarization induced in the atomic medium due to a linearly polarized probe field can be expressed as 

P - {i+x+Sp+ + e-X-Sp-) + c.c. , (5) 

where x± gives response of the medium corresponding to a± component of the electric field and, in general, is function 
of the fields £p±. In a previous work Agarwal et al |5[| studied the MOR under condition of saturation. However here 
we consider propagation of a weak probe and thus x± ^re independent of £p± . In the steady state the solution of the 
wave equation for Ep± gives the output field amplitude 

Eout {z^l) = e+Ep+ (0)e2"^>'>^+ + e_£p_ (0)e2-^'=-'^- . (6) 

For an x-polarized incident field say £in — xSq, 

£:p+(0)=£p_(0) = ^, (7) 

the output field in Eq.(6) takes the form 

Thus the output field Sout consists of both x and y-polarization components. The direction of light polarization 
rotates with respect to the polarization of the incident light. Experimentally the rotation is measured by measuring 
the intensity of transmission through a crossed polarizer at the output. In the above case the intensity of transmission 
through a y-polarized analyzer is given by 

J, ^ \{£out)y\'^ ^ 1 \2iiik^lx+ _ 1 2 (-g-> 

where the output intensity is normalized with respect to the input. It may be noted that for a resonant or near- 
resonant probe field, x± "^iH be complex and therefore due to large absorption by the medium, there will be large 
attenuation of output signal. However assuming that x± is real, i.e. absorption is small, we get from Eq.(8) 

-ii = tane* = tan7rfc;i?e(x- - X+). (10) 
For a commonly studied transition j — Q ^ j — \, the above susceptibilities are given by 

" vi^J {5±n-i-fy ^^^^ 

where 2hrL is the Zeeman splitting between the levels |j = 1, m = 1) and \i ^ \,m = —\)] 5 is detuning between the 
the probe frequency and the frequency ujiq of the transition j = 1 to the ground state j — Q (with zero magnetic field) 

(5 = wio-Wp. (12) 

In Eq.(ll), 27 is the decay rate of say \j = 1,711 — ±1) to the level \j = 0, m = 0). The quantity al gives the resonant 
absorption 

where n is the density of atoms and d is dipole matrix element for the transition. It is to be noticed that x+ = X- for 
zero magnetic field. Clearly to produce large MOR we have to make x± differ from each other to maximum possible 
extent. In next section we calculate the effect of a control laser on x±- 
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III. CALCULATION OF x± IN PRESENCE OF A CONTROL LASER 



We consider a model system [see Fig.l] involving say cascade of transitions \j = 0, m = 0) (level \g)) ^ [j = 1, m = 
±1) (level |1) and |2)) ^ \ j = 0, m = 0) (level |e)). This for example will be relevant for expressing ^°Ca. The 
probe Ep will act between the levels l^) and |1), |2). We assume in addition the interaction of a control laser Ec to 
be nearly resonant with the transition |e) ^ |2). For simplicity we drop the transition m = <-> m = 0. We thus 
assiimc the loss to m, = state by spontaneous emission could be pumped back by an incoherent pump. Let ujafj 
be the transition frequency between the levels \a) and The total Hamiltonian of the atom interacting with the 
control and probe fields is 

n = ?kOeg\e){e\ + ?lu;ig|l)(l| + hL02g\2){2\ 

-d-{Ec + Ep), (14) 

where both E^ and Ep are given by Eqs.(2) and (3). We make rotating wave approximation and thus we approximate 
the interaction part of the Hamiltonian by 

-d-{Ee + Ep)^ -Gi|e)(l|e-^"=*+''==^-G2|e)(2|e-^'^=*+''==^ (15) 
- 5i|l)(<7|e-''^''*+''=''^ - 52 1 2) (<7 16"''^''*+''=''^ + H.c, 

where 2Gs and 2gs represent Rabi frequencies of the control and the probe laser - 

Dei ■ £c _ dig ■ £p 

- - (16) 

In order to proceed further we will make following transformations on the off-diagonal elements of the density matrix 
- this transformation removes all explicit dependence on the optical temporal and spatial frequencies 

Pei = Peie-"^''+"'^', Pig = Pige-^-''*+^'=«'^ (17) 

In the following we drop tildes from the density matrix equation. However it should be borne in mind that the 
full density matrix in Schrodinger picture is to be obtained by using Eq.(17). On introducing the various rates of 
spontaneous emission, we can write the equations for density matrix elements as 

Pee = -2(ri + r2)pee + iGipie - zG*/9el + iG2p2e " «GSPe2, 
Pel = -(Fl + r2 + 71 + «(A - ri))pel + «Gi(pii - Pee) + iG2P21 ~ iglPeg, 
Pe2 = -(Fi + r2 + 72 + i(A + ^))pe2 + «G2(p22 - Pee) + iGipi2 - ig2Peg, 
Peg = -{Ti + r2 + i{A + 5))peg + iGiPig + iG2p2g - igipel - i92Pe2, 

pn = 2Tipee - 271P11 + iGlpei - iGipie + igiPgi - iglpig, (18) 

Pl2 = -(71 + 72 + 2iU)pi2 + iGlpe2 + igiPg2 - iG2Ple - iglPlg-, 

hg = -(71 + i{5 + ^))Pig + mifigg - Pn) + iG'lpeg - ig2Pi2, 

P22 = 2T2Pee - 272P22 + iG2Pe2 - iG2P2e + «fl2Pg2 - «52P2g, 
P2g = -(72 + i{5 - 0))/92g + ig2{Pgg - P22) + iG^Peg - igiP21, 
Pgg = 271P11 + 272P22 + WlPlg - WlPgl + ig2P2g " ig2Pg2, 

where the detunings are defined by 

A = (uJel — Ulc + ^) = (We2 — UJc — ^), W12 = 20., (19) 
5 = {wig - - f2) = {uJ2g —UJp + fl). 

InEq.(18) 2rj (27^) represents rate of spontaneous emission from |e) \i) (|z) \g)). However in further calculations 
we assume 71 = 72 = 7 for simplicity. The dipole matrix elements in (16) are given in terms of vectors (4) as 

Dei = -De+, De2=De-, 

dig = -de-, d2g = de+. (20) 
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Here D (d) denotes the reduced dipole matrix element corresponding to |e) <-> \i) <-!■ \g)) transitions. Clearly we 
also have 



9i = ^'92 

This should be kept in view to determine the component of circular polarization that connects various transitions. 

We next calculate the susceptibilities using solutions of Eq.(18). The induced polarization at frequency Wp will be 
obtained from off-diagonal elements pig,p2g- 

V = nTr{dp) = n{dgipig + dg2P2g + c.c.) (22) 

= n{dgipig + dg2P2g)e^''^''^^''"'-" + c.c. . 

The exponential factors in (22) come from the transformation in (17). As indicated earlier, pig and p2g will be 
computed to the lowest order in the probe field. Defining 

Pi9={^^^^^P2a^[^^s-, (23) 

and using (20) and (21), we can write (22) in the form 

V - {x+£p+e+ + x-£p-ne-''^''^''''' + c.c, (24) 

where 

We have been able to obtain analytical solutions for which we give below - 

+ _ 11 [\G2? + (7 + ^{5 - ^Wi +T2+ i{A + 5))] 



\G2\H1 + i{s + n)) + (7 + i{s - n)) [|Gi|2 + (7 + i{s + n)){Ti + ra + i{a + s))] ' 

h + (7 + as + n)){r, + r2 + i{A + s))] 

|Gi|2(7 + i{s - n)) + (7 + i{s + n)) [|G2|2 + (7 + 1(6 - n)){Ti + + z(a + s))] ■ 



(26) 



(27) 



These reduce to well known results in the absence of the control field Gi = G2 = 0. On substituting (26) and (27) in 
(9) we find the y-component of the transmitted field 



T -i 



al ,\ f al _ 

exp [ Jys-^ 1 -exp I Zys 



(28) 



As mentioned earlier, we choose the parameters in such a way that there is maximum asymmetry between x+ ^nd 
An important case occurs when, say, G2 = i.e. the control field is cr_ polarized {£c+ = 0, £c~ 7^ 0). Clearly x- 
or has the value in absence of the control field 

' = {i+Z-^)y ^^^^ 

where as x+ is strongly dependent on the strength and frequency of the control field 

+ ^ i-f{Ti+T2+i{A + 6)) 

^ |Gi|2 + (7 + i{5 + n)){Vi +T2 + z(A + 5)) ■ 



(30) 



For large |Gi|, both real and imaginary part of s+ will show Autlcr-Townes splitting. Note further that even when 
no magnetic field is applied = 0, 

s-{n = {])^ s+{^ = Q), for|Gi|7^0. (31) 

In this case we have control laser induced birefringence. A particularly attractive possibility is to consider the case 
when |Gi| is large and that we work in the regime of other parameters so that js+l ^ |s~|. Under such conditions 
we will have large MOR or large signal Ty. The experiments of Wielandy and Gaeta focus on the laser induced 
birefringence. 
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IV. NUMERICAL RESULTS ON LASER INDUCED BIREFRINGENCE AND ENHANCEMENT OF 

MAGNETO-OPTICAL ROTATION 

In this section we present numerical results to demonstrate how the MOR can be enhanced. Note that the parameter 
space is rather large and the results will depend on the choice of |Gi|, 1021, magnetic field, control laser detuning, 
and of course the probe laser detuning. We have carried out the numerical results for a large range of parameters and 
present some interesting results below. In all the numerical results we scale all frequencies with respect to Ti(= 7) 
and we take the absorption co-efficient al as 30. 

In Fig. 2 we consider the rotation of plane of polarization in the absence of the magnetic field. For no control laser, 
the rotation vanishes. For non-zero Gi, the medium becomes anisotropic (Eq.(31)), we obtain substantial rotation 
of the plane of polarization. The Fig. 2 also shows the real and imaginary part of x±- In Fig. 3, we show the results 
for non-zero value of magnetic field. We find definite enhancement in the magneto-optical rotation. The reason for 
this enhancement can be traced back to the flipping of the sign of Re x+ which is caused by the control field. In 
addition we can produce large rotation for probe frequencies in the neighborhood of the frequencies of the Autler- 
Townes components. In Fig. 3 we also show how the non-resonant control field can produce further enhancement. 
Our calculations also suggest some interesting domain in which the probe should be tuned to obtain large MOR. 
Application of the control laser permits us to obtain large MOR in totally different frequency regime. In Fig. 4 we 
show the changes in the transmitted signal if G2 ^ 0. The absorption and dispersion profile now exhibit a triplet 
structure. As noted earlier, intensity of the transmitted y-component of the signal depends on the asymmetry between 
X+ and When asymmetry goes down, as for example, in the third and fourth row (from (02) to (62)) of Fig. 4, 
then Ty can go down. In our model the control of x± is not independent of each other as we connect to a common 
final level. Clearly a large enhancement will be possible if x± could be independently manipulated. 

In conclusion we have shown how a control laser can produce birefringence as well as enhance the amount of 
magneto-optical rotation effect. Besides the control field can produce new frequency regions which show very significant 
magneto-optical rotation. 

GSA thanks S. Stenholm for pointing out Ref.[14] and for discussions. GSA is also grateful to Marian Scully for 
extensive discussions and insights on laser induced control of optical properties. 
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(b) 

FIG. 1. (a) The four-level model scheme (say of "^^Ca) having m-degenerate sub-levels as its intermediate states. The 
symbols in left hand side denote the energy levels of ''"Ca atom. 2Vi and arc the spontaneous decays, 2gi {2Gi) is Rabi 
frequency of the probe (control) field due to coupling of the intermediate state \i) with \g) (|e)). The detuning of probe (control) 
field from the center of |1) and |2) are represented by 5 (A). 2S7 is the Zceman split between the intermediate states, (b) A 
block diagram that shows the configuration under consideration. B defines the quantization axis z. The input probe Ein is 
a;-polarized and the control field is left circularly polarized. Both the fields propagate along z. After passing through the cell, 
output is observed through a j/-polarized analyzer. 
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FIG. 2. Plots to show the control field induced birefringence in absence of magnetic field (O = 0). Control field is left 
circularly polarized and is in resonance with \i) — > |e) transition (A = 0). In the above figure, first row refers to real and 
imaginary part s~ (Eq.(29)) which is independent of the O- control field. The frames (ttj), {bj) and (cj) represent Re s"*", Im 
s+ and Ty respectively. Three rows for j — 1, 2, 3 correspond to control field strengths Gi = 20, 50, 100 respectively. The 
X-axis represents the detuning of probe laser with respect to the position of the j = 1 level with zero magnetic field. 
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FIG. 3. Enhancement of MOR by use of a <t_ control laser with the intermediate levels being split by a magnetic field 
with Zccman splitting 2fl — 10. The graphs (aj) and (bj) represent the control field induced changes in Ro and Im 
respectively and (cj) represents the corresponding Ty, and j = 1, 2, 3 refer to the case of resonant control field (i.e. A = f2) 
with strengths Gi = 0, 20, 50 respectively. Clearly large fields result in significant enhancement and in addition new regions 
appear with large MOR. Further enhancement is observed by use of a detuned laser as seen in frames (04) and (cs) which are 
for Gi = 20 and control laser detuning A = —20, — 30 respectively. 
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FIG. 4. Plots represent MOR enhancement with an eUiptically polarized control (i.e. G2 7^ 0) and A = f2 = 5. Here we 

have considered the case of G2 = 10. The frames (flj), (&j), (cj), (dj) and (ej) represent Re , Im , Ro , Im a'^ and Ty 
respectively, and j = 1, 2, 3 correspond to Gi = 0, 20, 50 respectively. The off-resonant control field can also be advantageous 
(results not shown) as in Fig.3. 
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